The Use of an Electrical Potential Analyser 
for the Calculation of the Pressures on 
Lifting Surfaces 


Professor S. C. REDSHAW, D.Sc., Ph.D., M.I.C.E., F.R.Ae.S. 


(The University, Birmingham) 


Summary: The distribution of pressure over a lifting surface in subsonic potential 
flow is investigated by means of a three-dimensional electrical potential analyser. It is 
shown that by using the linear perturbation theory it is possible to obtain the pressure 
distribution for small angles of incidence, as well as the slope of the lift-incidence curve, 
easily and rapidly. A comparison with results obtained by vortex theory calculation 
shows a good agreement in spite of the relatively coarse network which was used for 
the experiments. There are no practical limitations to the use of a network having a 
finer mesh, and problems concerning curved and twisted surfaces may also be studied 

with facility by the use of the potential analyser. 


Introduction 


The distribution of pressure over a high aspect ratio lifting surface, immersed in 
a low subsonic stream of uniform velocity, may be calculated approximately by 
lifting-line theory. 


The determination of the pressure distribution over low aspect ratio swept and 
unswept plan forms presents considerable difficulty if analytical methods, based on 
lifting-surface theory, are used. 


An alternative method, which is described in this paper, is based on the use of 
an electrical potential analyser which is an electrical analogue machine in which the 
electrical potential, in a pure resistance network, represents the velocity potential of 
the disturbed flow in the vicinity of the lifting surface. From the measured potential 
the pressure can be calculated,- other aerodynamic properties being computed 
from the pressures in the normal manner. It is important to note that this method 
only provides a potential flow solution, although with circulation, for thin surfaces 
set at small angles of attack. 


Campbell” has discussed two electrical analogies for the determination of the 
pressure distribution over a lifting surface; the first analogy suggested by him is 
virtually the one used in this paper, but in the second analogy he suggests that the 
electrical potential should represent the acceleration potential of the fluid flow. In 
both cases Campbell proposed to apply the analogy by means of an electrolytic tank. 


Originally received June 1954. 
[The Aeronautical Quarterly, Vol. V, September 1954] 
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Malavard and Duquenne®, who used the velocity potential analogy, have 
recently studied lifting surfaces by means of an electrolytic tank; they preferred the 
velocity to the acceleration potential analogy except for certain special cases. It 
should be noted, however, that the acceleration potential analogy is valuable in 
cases where flaps are attached to the lifting surface and control hinge moments have 
to be calculated. 


A network will only provide an analogy to the finite difference form of the 
differential equation it is required to solve, as it is only possible to read the potentials 
at the discrete points provided by the nodal points of the net, whereas the electro- 
lytic tank provides a direct analogy to the continuous solution. However, the 
technique of using a resistance network is far simpler than that which has to be 
employed for an electrolytic tank. 


A detailed account of the evolution of the method, together with the recorded 
results will be found in a previous report™. 


Notation 


x,y,z Cartesian co-ordinates. Constant values of z define the plane of the 
tiers, z=0 refers to the plane of the master tier 


X,Y,Z Cartesian co-ordinates obtained by an affine transformation from 
co-ordinates x, y, Z 


U,V,W _ velocities along the x-, y- and z-axes respectively 
u,v,w perturbation velocities along the x-, y- and z-axes respectively 
® velocity potential 
@ perturbation potential, also electrical potential 
M Mach number 
p fluid pressure 
p, _ fluid pressure in the undisturbed stream 
Ap _ excess pressure defined by Ap=p — p, 
p fluid density 
q dynamic pressure 
C, pressure coefficient 
C, _silift coefficient 
angle of incidence 
lifting surface span 
lifting surface chord 


b 
c 
S lifting surface area 
A aspect ratio 
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2. Aerodynamic Theory 


2.1. Three-dimensional fields of flow 


If it is assumed that an incompressible fluid, in two-dimensional potential flow, 
possesses neither vorticity nor viscosity, then the flow may be described by means 
of a velocity potential function @ such that the velocity is equal to the gradient of ?. 
Owing to the continuity condition there also exists a stream function whose gradient 
gives the mass flow per unit length, normal to the flow. A velocity potential exists, 
whether the fluid is compressible or incompressible, but a stream function exists only 
when compressibility is negligible, irrespective of whether the motion is rotational 
or irrotational. 

In three-dimensional fields of flow the velocity potential still satisfies Laplace’s 
equation and is analogous to the electrical potential for the steady flow of electricity 
through a block of conducting material. The stream function, however, is no longer 
defined, except in the special case of axi-symmetrical flow. For this reason, when 
considering three-dimensional flelds of flow, the direct analogy must be used in 
which the velocity potential is represented by an electrical potential; the indirect 
analogy, in which the stream function is represented by an electrical potential, is 
inadmissible. 


2.2. Linear perturbation theory 


Within certain limitations, the linear perturbation theory considerably simplifies 
the analysis of some cases of potential flow. If a lifting surface is placed in a uniform 
stream which is moving with a velocity U parallel to the x-axis, say, the motion will 
be perturbed and the velocity potential will become 


P=Ux+¢(x, y, 2) : (1) 


so that » is the perturbation potential. 


The velocity at any point which was formerly (—U, 0, 0) now becomes 
(—U+u, v, w). The perturbation is said to be linear if u/U, v/U and w/U are 
small quantities of the first order, whose squares and products are negligible. No 
limitation is imposed on U, which may be either large or small, but the approxima- 
tion fails near a stagnation point where —-U+u=0, giving u/U=1. 


The assumptions implicit in the exercise of the linearised theory are 
(i) The lifting surface is represented by an infinitely thin plate. 

(ii) The camber of the surface must be small. 
(iii) Only small angles of incidence may be considered. 


(iv) The vortex sheet at the rear of the lifting surface remains in the same plane 
as the surface and runs immediately aft. 


These limiting assumptions define a lifting surface which would be suitable for 
high speeds and therefore, for this condition, the linearised theory will give 
good results. 
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The Prandtl-Glauert equation, satisfied by the perturbation potential is 


= 
ay? 4+ =() i (2) 


0*o 
— 2 
(1—M?) + 
where M is the Mach number of the undisturbed flow. 


As has been shown by Glauert“ and Prandtl®’, a distribution of potential 
satisfying Laplace’s equation will also satisfy the linearised compressible flow 
equation if the distribution 9 (x, y, z) is foreshortened along the direction of motion 
by the affine transformation 

Yxy 
Z=z2. 


Using (3), equation (2) transforms to Laplace’s equation 


ox? * oy? + 9: x (4) 


Hence, starting with a fictitious lifting surface longer in the x-direction than the 
true one, and calculating the potential distribution over this surface by incom- 
pressible flow methods, the correct dimensions and correct distributions of ¢ are 
obtained when the affine transformation is applied. 


The assumptions of the linearised theory allow the free stream component 
velocities V and W, in the y- and z-axes respectively, to be ignored in comparison 
with the velocity U along the x-axis. Now if p is the pressure at any point and p, 
is the pressure in the unperturbed stream, from Bernoulli’s theorem 


p+4p[(—U+u) +w*]=p, + @ 
Defining the excess pressure as Ap=p-p, ; ; : . (6) 


then, from (5), on the assumption that u, v and w are small, 


It has been seen that the neglect of the squares and products of the perturbation 
velocities and small angle of incidence leads to the result that u=U at a stagnation 
point. This entails infinite pressure at the leading edge except for the ideal angle of 
attack. Campbell has pointed out that there is no trouble in the absence of 
discontinuity in the streamline direction at the leading edge, for in this case (u=0) 
also the infinite velocity may be removed, except at the stagnation point, by round- 
ing the leading edge. At a finite lift the stagnation point moves to the lower surface, 
but if the thickness, curvature and angle of incidence are diminished, the domain in 
which the error is considerable shortens to the immediate locality of the 
leading edge‘. 
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The Joukowski condition, that the air speed at the trailing edge must be finite, 
has to be satisfied in order that the flow of an ideal fluid may approximate to that of 
a real fluid. This postulates that the streamlines at the trailing edge of a thin 
surface must be continuous in direction with no infinite velocities. 


The incidence of the lifting surface is given by the ratio of the vertical velocity, 
in a region not affected by the wing, and the horizontal velocity, thus 


W 
— Ap _ 2u 
and from equations (7) and (8), (9) 


As the pressure on the lower surface of the wing is equal and opposite to the 
pressure on the corresponding point of the upper surface, the lift coefficient is 
given by 


b/2 
=4|4 2 Pas (0) 
0 
b/2 c 
0a da 


3. The Electrical Analogy 
3.1. The theory of the analogy 


If fluid flow is identified with electrical flow, there is a direct analogy in which 
the velocity potential is represented by an electrical potential. Now as 


u= —09/0x . ‘ & 
and, at a point remote from the surface, 

W = — 00/0z ‘ ‘ (13) 
then equation (9) may be rewritten as 


Ap 2/0x 
qu 


(14) 


where » denotes either velocity or electrical potential. Then from equations (11), 
(12), (13) and (14) 
b/2 


—b/2 


oc, 4 
Sdp/dz 


since » is zero at the leading edge. 
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Fig. 1. 


3.2. The three-dimensional potential analyser 


A description of the three-dimensional potential analyser has already been 
given. The analyser consists essentially of a cubical mesh of resistors arranged in 
nine main tiers with one auxiliary tier, each tier consisting of a square array of 
25 x 25 sockets interconnected with resistors to form a square mesh. The nodal 
points on each tier are connected to the adjacent tiers through inter-tier resistors 
having the same resistance value as the tier resistors. The individual mesh 
resistances consist of special precision woven resistors having a tolerance of +1 per 
cent. on a nominal resistance of 200 ohms. 


Resistance elements were not assembled on the master tier, but post office type 
terminals were used instead, these terminals being connected to the inter-tier 
resistors in the normal manner. Resistors, depending on the nature of the particular 
experiment in hand, were attached to the terminals of the master tier as required. 
The resistors on two adjacent edges of each tier were doubled in value and the 
inter-tier resistors connecting these edges were also doubled, while the inter-tier 
corner points enclosed by these edges were quadrupled in value, thus providing 
“selvedges” for problems involving a field of single or double symmetry. The 
reason for using the doubled resistors is that the net can be considered to be 
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split at the plane of symmetry. When resistance elements were added to the master 
tier, which is a plane of symmetry, they were also doubled in value. A view of the 
instrument is given in Fig. 1. 


3.3. The realisation of the analogy 


The potential analyser can then be used to realise the analogy by 


(i) Slitting the bottom tier in its own plane over an area corresponding in 
shape to the aerofoil considered, 


(ii) Short-circuiting the area surrounding the model and short-circuiting the 
whole of the top tier, and 


(iii) Applying a known electrical potential between the short-circuited top tier 
and the short-circuited margin of the bottom tier. 


In practice the bottom tier is not actually slit because, as the flow is applied 
normal to the model, this tier may be considered as a plane of symmetry. It is, 
therefore, only necessary to double the value of the resistances over the area which 
represents the model. This analogy is a three-dimensional form of the direct 
analogy of Taylor and Sharman’, When flow with circulation is being investigated, 
the Joukowski condition at the trailing edge has to be satisfied; this may be realised 
by ensuring that points immediately behind the trailing edge shall be at the same 
potential as corresponding points along the trailing edge. On the potential analyser 
this is effected by short-circuiting lines along the x-axis aft of the trailing edge, and 
raising the potential of each individual line so that its potential is the same as at 
the corresponding point on the trailing edge, one mesh separation away from it. 
Thus in the wake, as ¢ is independent of x, 


Ww = TR 


OW 

CX 

where suffixes “w” and “TE” denote conditions in the wake and at the trailing 

edge, respectively. 


The second of these conditions applies because the pressure, which is propor- 
tional to 0¢/0x, vanishes on both sides of the wake and is, of course, zero at the 
trailing edge. must be continuous at the trailing edge, as any discontinuity would 
cause an infinite velocity there. 


4. Experimental Procedure 


The analogy, described in the previous section, was applied in the following 
manner; the plan form of the lifting surface under consideration was set up on the 
master tier by attaching resistors to the terminals over the area concerned, the rest of 
the tier being short-circuited with the exception of the area defined by the wake. 
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Fig. 2. 
Electrical circuit of analogy for flow with circulation. 


edge, were shorted to the boundary of the tier. Thus the area from one unit behind 
the trailing edge to the boundary was shorted by strips running in one direction only; 
each shorting strip was connected through a rheostat to a potential divider. A 
diagram illustrating the circuit used is shown in Fig. 2. In each experiment only 
half the lifting surface was set up on the master tier, the axis of symmetry of the 
surface being arranged to lie along one of the selvedges of the tier. To avoid 
channel restraint, due to the proximity of the boundary, the span of the surface 
occupied about half the tier width. 


The experimental procedure was to apply a unit voltage between the short- 
circuited points of the auxiliary tier and the short-circuited area on the master tier 
and then, by means of the potential divider and rheostats, to adjust the potentials 
on the shorting strips so that they agreed with the potentials at corresponding 
points on the trailing edge of the surface, thus satisfying the condition for circulation 
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that the potential at a point immediately behind the surface shall be at the same 
potential as the corresponding point on the surface trailing edge. The adjustment 
of the potential on any strip upset the values of the potentials on the neighbouring 
bars, and the setting-up procedure developed into what might be termed an electrical 
iteration process. However, with a little practice, it was found that the setting up of 
an experiment was quite rapid. When the set-up had been completed the potentials 
on the lifting surface, as defined on the master tier, were measured in the normal 
manner by the use of a Muirhead Type D.72 A potentiometer. All readings were 
taken to the nearest millivolt, the limit of accuracy of the potentiometer. It should 
be noted that the only potentials which have to be measured are those at the nodes of 
the lifting surface on the master tier. 


5. Results 


The following three lifting surfaces were investigated : — 

(i) Rectangular, constant chord, aspect ratio 6. 

(ii) Delta, 45° sweep at leading edge, aspect ratio 4. 
(iii) Swept-back, 45° sweep at leading edge, constant chord, aspect ratio 3. 


The electrical potentials, directly analogous to the velocity potentials, which were 
determined directly from the experiments were plotted for chordwise sections of 
the surface; these curves were graphically differentiated by means of an optical 
differentiating device to give values of 0¢/dx. Then, from equation (14), values of 
Ap/(q2) were computed; these results are represented pictorially in Figs. 3, 4 and 5 
and plotted in coefficient form in Figs. 6, 7 and 8, where they are compared with the 
results of vortex calculations®”’ supplied by Dr. D. Kiichemann of the Royal Aircraft 
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Fig. 4. 


Establishment. The experimental results agree well with the calculations except in 
the region of the leading edge, where theory predicts an infinite pressure which 
cannot be obtained in actual fluid or electrical flow conditions. 


Values of 0C,,/0a for each surface were calculated from equation (15) and 
compared with theoretical values supplied by Dr. W. P. Jones of the National 
Physical Laboratory. This comparison is given in Table I. 


The calculation of 0C;,/0a from the experimental results was extremely rapid; 
it was only necessary to measure the potentials at the trailing edge, plot the spanwise 
distribution from these figures, integrate the curve and evaluate equation (15). When 
making these calculations, however, it was difficult to assess the true area of the 
aerofoil. The difficulty arose because the net used was too coarse for the true 


TABLE I 
Lifting Surface ‘Theoretical xperimental- 
Rectangular (A = 6) 4:26 4-80 
Delta (A=4) 3°47 3°22 
Swept-back (A =3) 2:75 2°75 
172 


: 
< > 
— 
— 
1 


ELECTRICAL POTENTIAL ANALYSER 


Fig. 5. 


boundary position to be precisely located. The actual locus of the boundary must 
lie between lines joining the nodal points and the assumed boundary, which was 
short-circuited on the surface set up on the analyser, and lines joining the nodal 
points one mesh distance within the surface from the short-circuited boundary. For 
the purpose of calculating 0C,,/02, which involves wing area, the true boundary 
was assumed to lie along lines one half mesh distance inside the model from the 
short-circuited boundary. This difficulty would have been obviated if the mesh 
separation of the tier had been finer. 


6. Conclusions 


The potential analyser is a particularly easy and rapid instrument to use, even 
by unskilled personnel. The results of the experiments are in good agreement with 
theory, in spite of the disadvantage of having a relatively coarse net. There are no 
practical limitations to the building of an analyser having a much larger master 
tier, with a smaller mesh, and a considerable constructional simplification would be 
possible, for the type of problem considered in this report, as only the master tier 
need be accessible for scanning. It is important to note that curved and twisted 
surfaces may also be studied by the use of the potential analyser without difficulty. 
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A Note on the Numerical Solution of 
Fourth Order Differential Equations 


L. C. WOODS, D.Sc., A.F.R.Ae.S. 
(Department of Applied Mathematics, University of Sydney) 


Summary: An old numerical method of solving fourth order differential equations 
is put in relaxation form. The higher order correction terms are included and the 
technique is illustrated by an example. The method has the advantage of being more 
rapidly convergent than the usual relaxation procedure for fourth order equations. 
Some comments are made on the numerical solution of the viscous flow equation. 


1. Introduction 


Southwell’, Fox" and Thom? have each described methods of solving two- 
dimensional fourth order differential equations of the form 


of o4 
where ox! +2 Ox*dy? + 


The first two authors used “relaxation” while Thom used an iterative process he 
termed “ squaring.” There was one other important difference in the methods. This 
was that Thom used an intermediate function ¢ (x, v), such that 


V*w=(, and hence V*(=k(x,y), . (2) 
and solved these two equations of Poisson’s form instead of the single equation (1). 


In this paper a relaxation form for Thom’s method is developed, and details are 
given of the treatment of the boundary conditions most commonly experienced. The 
higher order correction terms are also given, enabling a greater accuracy to be 
achieved on a mesh of a given size. A simple example of equation (1) is solved to 
illustrate the technique. One-dimensional fourth order equations can be solved 
similarly. 


The outstanding advantage of solving the simultaneous equations (2) instead of 
equation (1) lies in the very much greater rate of convergence of the relaxation 
process for Poisson’s equation, compared with equation (1). Fox‘ has considered 
this point in some detail. 


Originally received July 1953. 
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Notation 


h mesh size 
a difference operator corresponding to the Laplacian V? 
k the right hand side of Poisson’s equation 
the operator 0°/0x dy 
M_ modulus of a conformal transformation 
R”,R’ relaxation residuals 


0/0n normal derivative on a boundary. 


2. Equations for the Residuals 


If f,, i=0, 1,2, . . . are the values of a function f (x, y) at the mesh points shown 
in Fig. 1, then it can be shown that“? 


4 6 
where /: is the mesh size, 
a 
and Af,= fi-4f,. 
t=1 
If an auxiliary function wv is defined by 
so that from (1), h?V4w= V*v=h’k, (5) 


then, with the aid of (4) and (5), an application of (3) to the function w results in 


Similarly from (3), (4) and (5) 


Av, =h*k, + 12 (h?V?k, —2D‘v,) + 360 V'k, +O (h’). (7) 


An application of the operator A to equation (6) yields 


A’w, + {(h?V°*k, Ak,) 2D* (vy + Aw,)} O (h*), 
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and so, from (3) applied to k and w, it follows that 
2 2 2 4 8 
From Fig. | it is found that i 
A 8 12 


Equation (8) can be established directly’, but the value of the derivation given 
here is that it shows that the correction term h* (V°k, —2D*v, /h?), is still obtained 
even if the O(h‘) terms are ignored in (6) and (7). If equations (6) and (7) are 
solved by the usual relaxation method: for Poisson’s equation, taking into account 
the fourth order terms, then the sixth order term in equation (8) will be correctly 
allowed for. 


It is thus necessary to solve the simultaneous equations 


~ 


Aw, 2D'w,), | 


2D). 


and Av, =h'k, 4 12 


A method of solving single equations of this type has been described in detail in 
other publications”:®’. As a first approximation the equations 


Aw,=%,, Av, =h'k,, . » 


are solved throughout the field by the usual relaxation technique’. Residuals R" 
and R° are calculated at point 0 by 


Aw, — Vos 


R,* = Av, 
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From (12) and similar equations for R;, i=1, 2, 3, 4, it follows that 


OR: = i= 3, 4; OR, wih, 

Ow, Ow, 

oR v oR v oR w eo % (13) 


Fig. 2 shows a suitable method of recording the values of w and v and their residuals 
at a typical mesh point. When residuals defined by (12) have been sufficiently 
relaxed by (13), the residuals defined by 


h‘ 


R,”= Aw, 6 D‘w, 12 
(14) 
R,”= Av, — h*k, + 
0 6 12 0° 
are computed and relaxed. If an operator A’ is defined by 
A’w,= wi- 
t=5 
then (h*/6)D‘w, can be written? 
h* 
— D‘w,= wi-2 wit+4w,} +0 (A) 
6 Stns i=1 
™ (A’w, —-2Aw,) +0 (h'). 
h* 
Similarly = D‘v,= (A’v, —2Av,) +0 (A), 
and so (14) can be written without loss in accuracy in the form 
R,”= (4Aw, + A’w,)— k 
an 6 0 0 “oe 12 0° 
(15) 
R,’= (4Av + A’v,)—h*k, 
“0 “0 1 0 12 


It is convenient, but not essential, to relax these residuals by the simple relaxa- 
tion patterns (13), but of course this introduces small errors, which are eliminated by 
recalculating the residuals from (15) after each complete relaxation using (13). The 
principle involved here occurs commonly in relaxation—the residuals are calculated 
accurately, but simple and approximate relaxation patterns are used. 


The method of Fox and Southwell’) is to solve the single set of difference 
equations given by (8) and (9), but ignoring the correction term. In a later paper’? 
Fox considered the correction term in a different form. 
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3. Boundary Conditions 


Boundary conditions of common occurrence are values of w and dw/dn (the 
normal derivative) specified on boundaries which in general are curved. Only these 
boundary conditions are specifically treated in this paper, but the method can, 
if necessary, be applied to other types of boundary condition. From the computor’s 
point of view the simplest boundaries are those composed of straight lines parallel 
to the x- or y-axis and passing through the mesh points. Curved boundaries can be 
reduced to these simple boundaries in two ways, namely, (a) by extrapolating the 
boundary conditions to external mesh points (as the calculation proceeds) to obtain 
pseudo-boundary conditions on boundaries parallel to the x- or y-axis (see Ref. 1), 
or (b) by a conformal transformation of the xy-plane to eliminate the curved 
boundaries (see Ref. 8). The former method has—at least for viscous flow problems 
—disadvantages compared with the conformal transformation, the theory of which 
now follows. 


It is required to solve 
when w and ow/0n are given on curved boundaries. The first step is to make a 
conformal transformation to an 28-plane so that the boundaries become straight and 


parallel to either <=0 or B=0. 


Let M be the modulus of this transformation, so that 


+ (sy. 
Then (16) becomes w) =k(a,B). . . (17) 


If an auxiliary function v is now introduced, defined by 


Vv 
2 
V7. (18) 
where / is the mesh size in the 28-plane, then from (17) 


Corresponding to equations (10) it is found that 
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Now consider a boundary AB parallel to 8=0 (Fig. 3). On AB, 


=s (a), J ‘ (20) 


w=r(q), 
say, where r and s are known. Now (see Fig. 3), 


war th( Se) + + (21) 


where, from (20) r,=W,, 5,=(0w/0f),, and from (18) and (20), 
Ga),* 


where some use has been made of (19). These results enable (21) to be written 


Vo h(a v v 
(22) 
h? ( d’r h® ( (dr h*k, 
where E,=r,ths,— % x2), 24M,*’ 
is a known function. Also 


When {(0/08)(v/M7?)}, is eliminated from (22) and (23), it is found with the aid 
of (19) that 


which gives a comparatively accurate boundary condition for v. If second order 
terms are ignored, (24) reduces to 
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4. Some Comments on the Viscous Flow Equation 


In the ~8-plane the equations of viscous flow take the form? 


1 (ay a ay 
2 3 > 


where v is the kinematic viscosity, ¥ is the stream function and (¢ is the vorticity. The 
boundary conditions are 


(28) 


Equations (26) and (27) can be combined into one fourth order equation for y. 
but of such complexity that its direct numerical solution would be almost impossible. 
After lengthy calculations, Thom managed to solve the simultaneous equations (26) 
and (27) by his “squaring” method for several low values of the Reynolds number 
for the flow about a circular cylinder’. The recent development of high-speed 
computing machinery has given this early work of Thom a new interest. 


The relaxation process described in Sections 2 and 3 can be applied to this 
problem. It is of some interest to write down the boundary equation for (, since 
it differs from the one given by Thom. The value of k (see equation (16)) for 
viscous flow is plainly 


1 (a4. _ 
v\0x dy dy dx/’ 


which vanishes on the boundaries. Thus for this example E, in (25) reduces to y¥,. 
Thus when w and v are replaced by ¥ and h’¢ respectively, equation (25) becomes 


2 2 


which should be used in place of equation (8) of Ref. 8. 


A final point that should be noted is that, near the boundaries, (27) is the 
dominant equation of the motion, whereas away from the boundaries (26) dominates. 
In terms of residuals this means that RY and R‘ have a relative importance which 
varies with distance from the boundary. If the usual rule “eliminate the residuals 
in order of size” is to be applied, then R‘ must first be multiplied by a weighting 
factor proportional to the value of ¢ at the point in question. The author believes 
that if the residuals are not weighted in this way the relaxation process may become 
divergent in certain cases. At large Reynolds numbers there is another possible 
cause of instability in the relaxation process, which was implied in Ref. 8. It is that 
the flow may be basically unsteady, and the assumed differential equations will not 
adequately describe it. This may show up in the relaxation process as divergence 
or instability. 
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5. An Example 


The simplest example of equation (1) is probably that of finding the deflection of 
a thin square plate, clamped at its edges and loaded uniformly on one face. It has 
been calculated previously by Fox and Southwell” and later with greater accuracy 
by Fox". 


The boundary conditions are w=0, 0w/én=0 on x=+1, y= +1, and the 
symmetry makes it necessary to consider only one octant of the square. The value of 
k in equation (1) will be taken to be 25,600 so that the results may be compared 
easily with those given by Fox. In this example equation (24) becomes 


The solution in the cases h=4 and } is set out in Fig. 4. The values closest to the 
horizontal lines (see Fig. 2) are those obtained by solving equation (11) and using 
the comparatively accurate boundary conditions (30), while the outer figures are 
obtained by solving the more accurate equations (10), again using equation (30) on 
the boundaries. In the case h=4, Fox" obtained the figures given in circles to the 
right of the mesh points by using a method due to Hencky referred to in his paper. 
Fox’s relaxation solution (obtained by relaxing the biharmonic equation directly, 
and not shown in Fig. 4), the author’s solution, and the results from Hencky’s 
method differ at most by one unit. 
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0 
0 
i 
| i -226 26 0 
| 
‘ | 33 77 
: 177] 0 
: 188 0 188 70 (9) 
| -151 840 
=56 33 208 
| 852 | 37 309 
414 278 102 0 
389 4 | 405 274 103 
| | 178 “104 42 297 
| 
Gis) 463 310 114 
h = -507 Gis) 453 305 G1) 14 @ Q 
-214 -189 -106 54 327 
-226 -200 =ttS 49 330 
= 14 
Fig. 4. 


The labour involved in obtaining the solution shown in Fig. 4 is perhaps a 
little greater than that involved in solving Poisson’s equations with double the 
number of mesh points and given boundary values. It appears from Part B of Fox’s 
paper that this is very much less labour than that involved in solving the single fourth 
order equation directly, on this mesh. 
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A Theory of Uniform Supersonic Flow 
Past a Thin Oscillating Aerofoil at 
Appreciable Incidence to the Main Stream 


GEOFFREY L. SEWELL, M.A., Ph.D. 
(University College, Hull)* 


Summary: A number of authors (e.g., Refs. 1 and 2) have dealt with the problem 
of uniform two-dimensional supersonic flow past a thin oscillating aerofoil, whose angle 
of incidence to the main stream is small. In this special case, the effects of vorticity are 
negligible, while the fields of flow on either side of the aerofoil may be treated on the 
same footing. 

The object of this paper is to deal with the more general case in which the aerofoil 
performs small oscillations about a fixed position, in which its incidence to the main 
stream is appreciable. The oscillations modify both the Prandtl-Meyer expansion on 
one side of the aerofoil and the shock compression on the other side of it. Both these 
effects are treated by a linearised theory. It is found that vorticity is generated down- 
stream of the shock line, and full account is taken of this. 

It was shown in Ref. 2 that the flutter derivatives for a slender aerofoil are the 
same as those for a flat aerofoil of the same chord, oscillating in the same manner, and 
therefore only the case of the flat aerofoil is treated in this paper. 

The method is sufficiently general to cope with all types of flutter, although solutions 
in closed form are actually obtained only for the cases when the frequency of oscillation 
is small. 


Notation 


OA mean position of the chord of the aerofoil, O being the 
nose and A the tail 


BO line of direction of impinging main stream — 


OS _ shock line on the upper side of the aerofoil when it is 
stationary 


OM,,OM. Mach lines through O on the lower side of the aerofoil 
between which the Prandtl-Meyer expansion takes place in 
the stationary case 


R,, R,, R., R,; regions between OS and OA, BO and OM,, OM, and OM.,, 
OM, and OA, respectively 


Oxyz axes of reference for flow in R,, Ox being OA and Oy lying 
in the plane of flow 


k unit vector in direction of Oz 


*At present at Marischal College, Aberdeen. 
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Fig. 1. 


Ma 


V (magnitude V), po’, Mo’ 


u, (magnitude U4), Po, Pos Qos 
u, (magnitude u,), p,, p,,a, and p, 


u(=(u, v, 0)), p, p 


n,, 


angle between BO and OS 
angle between BO and OA 
angle between BO and OS 
length of OA 


shock curve on the upper side of the aero- 
foil when it is oscillating 


velocity, pressure, density, speed of sound 
and Mach angle, respectively, in the 
impinging uniform stream 


corresponding quantities in R, when the 
aerofoil is stationary 


corresponding quantities in R, when the 
aerofoil is stationary 


perturbation velocity, pressure and den- 
sity, respectively, in R,, when the aerofoil 
is oscillating 


point on = 
foot of perpendicular from P to OS 


unit vectors along the normals to OS (at 
any point) and = (at P), respectively, their 
directions being into R, 
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t, unit vector along OS 
€ = P,P, in direction of n, 
s = OP, 

Os Ox dy 

dy 


m,’) normal and tangential components, re- 
spectively, of the velocity of flow relative 
to OS immediately upstream of it 


(1,, m,) normal and tangential components, re- 
spectively, of the velocity of flow relative 
to = at P, immediately upstream of it 

vy frequency of flutter 
v, (=ae™), 6(=be™), w(=6=ivbe™) flutter velocity component of the nose 
perpendicular to OA, angular displace- 
ment and angular velocity of the aerofoil, 
respectively, a and b being constants 


F,G __ resultant and pitching moment about O, 
respectively, of the forces on the upper 
side of the aerofoil due to its oscillation 


1. Modified Prandtl-Meyer Expansion on Lower Side of 
Aerofoil 

On considering any point of the aerofoil whose mean position is A, (on OA), it 
is seen that it sets up waves, the effect of which is nil upstream of the line A,A, 
(parallel to OM,). Hence the effect of the oscillation must be zero in R, and R,. 
The flow in R,, moreover, is seen to be that of a uniform stream of velocity u,, 
pressure p,, and density p,, impinging on a flat aerofoil which is oscillating about 
OA; and this solution has already been obtained (see, for example, Ref. 1). 


A 


d 
e 
e 
A, 
R 
R, M, 
\ 
M, 
Fig. 2. 
187 


GEOFFREY L. SEWELL 


2. The Flow in Ro 


The problem of the flow in R, is now treated by perturbation theory. The 
linearised equations of flow in this region are 


+pdivu=0,. . . . . 

grad p+ py (= on) 70, ‘ (2) 
and (< +U, =) (p-a,?p)=0. (3) 
Hence, by equations (1) and (3), 

+pazdivu=0; . . . . 


while, by equation (2), p and u may be expressed in the forms 


oH oH 
(A) 
and u=grad H + curl (ck), 
On substituting these formulae for u and p in equation (4) it is seen that 
Moreover, since v is the frequency of oscillation, 2 and H may be written as 
x 
a=a(y)expiv =), « « 
and H=$¢(x, y)expiv (1 - = sec? ‘ (8) 
where (by equation (6)), 
0? vy? 
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The shock conditions (see (2)) are 


u+u,- V= 


, 2 2 
and P+Po- Po = ay *) on OS, 


(B) 


!, being a variable because of the motion and curvature of =. The velocity of the 
shock at P is (0€/dt)n, and therefore 


L=(v- n= v.n— sin (o,- 26) 


=V sin 0, (Vcosa, + - (V cos, + 


while, since (6, — 0£/0s) is the inclination of the tangent to = at P to the x-axis, it is 
easily seen that n=n,—(0£/0s)t,. Hence the shock conditions (B) may be 
rewritten as 


2 0€ 2 ) 


and thence (by equations (A) and the shock conditions V cos 4, =u,coso and 
Po ly’ = poll, Sino) 


OH Oa _ _2 
ds On ds 
and + u, sino (u,cos + 37) OS. 


The boundary condition at the aerofoil (see Ref. 2) is 


v=(a+bu,+ivbxye“ ony=0. (10) 
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It will be convenient to arrange that <=0 on y=0, and to do this we observe that, 
by equations (A) and (7), 


u=grad H + curl { 


=grad [ - i2(O)exp { - ] + 
0 0 


+ curl exp iv (1 ~ *) k| 


Now -ia(0)exp [ +ivt| 


is a solution of equation (6), and so may be incorporated with H, while 2 (y) may be 
replaced by «(y)—2(O0)e~'”/“» which is equivalent to making (0) zero. In view 
of this, equation (10) may be rewritten as 


on y=0,.°. =. . (iD 
or, by equation (8), 
=(a+ bu, + ivbx) exp (= ) on y=0. 


By eliminating « and € from (C), a shock equation may also be derived for 
H(or ¢) only. To do this, 02/0n is first expressed in terms of 0a/0s and a. 


Since (by equation (5)), Ui, + 
da ( ag 26) 


Hence, after substituting this expression for 02/@n in the first equation of (C) and 
then eliminating € and from the resultant equations, it is found that 


A,H=0 on OS, ‘ (14) 
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where 
[ \( Os =) Ox =) 
1,7? +1 u,cose + +uU, Sine Uy + + 


and, with the aid of equations (8) and (14), 


where A is a differential operator of the third order in x and y only. 


The pressure distribution on the aerofoil can now be obtained by Riemann’s 
method. @® is a function satisfying the equation L?=0 (by equations (9) and 


(15)) and the boundary condition OS. If now we put V=J, /u,) sec? ) 


where (?=(x x,)? sin? », y,)? cos? then LYW=0 also. Hence, if P, Yo) 
lies on OS, and if PJ, PK are the Mach lines through P when tne aerofoil is station- 
ary, J and K lying on Ox, and if A is the area enclosed by the triangle PJK, then, by 
Green’s Theorem, 


= | [ (WL® dxdy=0, 


7 oO J K A 
| Fig. 3. 
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K 
ov 
which reduces to (o - dx=cot u,(Py+Px), (16) 
J 
X, Sin (4) — 7) ) 
J being the point 0 


( X, sin (u, +o) 0 ). 


and K the point sini 


Equation (16) is one involving the values of » and its derivatives up to the 
fourth order on y=0. 09/0y is known on this line (equation (12)), and therefore 


ap Op ) 
and hence By? and since L By =0 


are also determinable on this line; while, if ¢=f(x) on y=0, we can determine 


2. 2 

Ox’ Ox?’ dy?’ axdy?”’ dx*’ ax?dy’ 
(with the aid of equation (9)) in terms of f(x) and its derivatives, on y=0. This 
means that equation (16) may be transformed to an integro-differential equation for 
f (x), from which the pressure distribution on the upper side of the aerofoil may be 


arbitrary constants of integration. 


determined, although equations (C) will have to be invoked in order to evaluate the 


3. Solutions when V is Small 


When v is so small that its square may be neglected, equation (9) reduces to 


and equation (12) to 
=(a+ bu,)+iv [ sect xon y=0... (18) 
0 


This suggests a solution for ¢ which is composed of linear and quadratic terms 
only, the latter being of the order of magnitude of v; and since x, y, xy and 
x? +? cot? », are solutions of equation (17), we put 


¢=qx+(a+ bu,) y+iv [o+ {e+ 
0 


xy + ivr (x? + y? cot? 


where qg and r are unknown constants, this equation satisfying equations (17) ana 
(18). Hence (by equation (9)), 


[ ax+ (a+ bu,) ivbxy + ivr (x? + y? cot? - qx? sec? (19) 
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Further, we put E=(€,s+ivé,s?)e™ . ‘ 


and a=u, Sino on OS, . (21) 
2,, @, €, and €, being unknown constants. Hence, by equation (9), 


= = cosa on OS. 


Substituting these formulae in equations (C) and equating coefficients, gives six 
simultaneous linear equations for g, r, ,, 2, €, and €,, from which 


q=-u bu,)sin 2x (1+ 9% 008? o+ sino}. 


and 


a+ 


r(K coso +2u, sin? cos cot? + sin + 2u, sin cos? + 


+(a+bu,) sino [2 sec? cos” + sin? — 2u,K~-' cos? o (a; 1,7 + cos? =0, 


where K= (1+$ Me cos? + sin o. 


Thus we may say that 


q=4,a+ 
(D) 


r=r,at+r.b, 
where q,, q2, r, and r, are known constants. 


= + = , so that, on y=0, 


P= ~ Po\ Ot 
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Hence, 


F= PdX= — Po (1 — 2 sec? — 
0 


5 Pe (2ratly +4 (1 ~2 sec* c*6 


[2r,u, (1 -2sec? p,)] 


1 
and G= | xpdx= Poo Vo 6 Po 
0 


5 ~ py + 4a (I 28ec* 


from which the flutter derivatives may be obtained. 
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Flame Stabilisation in High Velocity Gas 
Streams and the Effect of Heat Losses at 
Low Pressures 


D. B. SPALDING, M.A., Ph.D. and B. S. TALL, B.Sc., B.E.(Sydney) 
(Engineering Department, University of Cambridge) 


Summary: The first object of this paper is to show the concordance between three 
different theoretical approaches to flame stabilisation and to compare published experi- 
mental data with the theoretical predictions. Good agreement between the theory and 
experimental results is obtained. The second object is to show that, at low pressures, 
the predictions of these theories must be qualified in view of the increasing importance 
of heat losses from the flame region. It is shown that there must exist a minimum 
pressure below which a flame cannot be maintained in a given combustible mixture, 
irrespective of the rate of flow of the mixture, unless heat losses are entirely eliminated. 


An attempt is made to calculate the value of this low pressure limit for the fuel-air 
mixtures used in aircraft engines, particular attention being paid to heat losses due to 
thermal radiation from the burned gases. Although the basic data of the calculations 
are somewhat uncertain, and some of the assumptions may be open to question, the 
effect appears to be of practical importance to engine designers. Some light is also 

thrown upon the existence of inflammability limits. 


1. Introduction 


The growing importance of the gas turbine and jet engine during recent years 
has focused considerable attention upon the mechanisms of flame stabilisation in 
a gas stream, and a number of semi-quantitative theories have been put 
forward":?-*-4.5), Three distinct approaches to this problem are considered here, 
and it is shown that, under corresponding conditions, they each yield the same result. 


Notation 
A _ entering gas mass flow rate 


a___ surface area of the primary reaction zone through which gases enter 
this zone 


B’ _ collision factor 

b exponent of d in Vgo Hd” 

c specific heat 

d__ characteristic dimension of the stabiliser 
E activation energy of a reaction 


H calorific value of the fuel 
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.. constants 


thermal conductivity 

molecular weight 

fractional mass concentration of fuel in a mixture 
order of the chemical reaction 

static pressure 

heat losses per unit mass of gas mixture 

heat release rate per unit volume 

universal gas constant 

radiation heat loss factor 


laminar adiabatic flame speed of the gas mixture referred to the 
temperature of the unburned state 


temperature 

average normal velocity of the gases entering the reaction zone 
velocity of gas stream at blow-out condition 

volume of reaction zone 


equivalence ratio of the gas mixture for lean mixtures; equal to 1-0 for 
rich mixtures 


reaction rate factor 

Peclet number based upon stream velocity 

rectangular space co-ordinates 

velocities in the x, y, z directions 

air-fuel ratio of a combustible mixture by weight 

thermal diffusivity 

constant 

mass fraction of the fuel burned in the primary reaction zone 
density of the mixture - 

represent functions of a variable 


Suffixes “b” and “u” denote the fully-burned and unburned states. 


2. Flame Stabilisation Theories 


2.1. Homogeneous reaction zone 


Longwell) considered the eddy region in the wake of a bluff body stabiliser and 
assumed it to be a region of intense mixing and so to be homogeneous in gas 
composition and temperature. By equating the difference between the masses of 


fuel entering and leaving this zone to the mass of fuel consumed in the chemical 
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reaction therein, he showed that for a bi-molecular reaction, equilibrium conditions 
are characterised by an equation of the form 


(1 —c)(1—ey) { 1 1) -2 
vp? R°T? M(a/f) 295 
where v= volume of the reaction zone (v replaces Longwell’s V which will be 


reserved for velocity), 


T = adiabatic flame temperature resulting from burning fraction < of the 
fuel, 


B’=collision factor (assumed by Longwell to be independent of 
temperature), 


and the other symbols are as given in the Notation. 


By the use of this equation Longwell showed that, at the condition correspond- 
ing to high-velocity extinction of the flame, the right hand side of equation (1) is a 
maximum, and any increase in the function A/(vp*) above this maximum leads to 
sudden extinction of the flame. Since the heat release rate per unit volume of the 
primary reaction zone is equal to (AH /v) (m;,—m;), where m,, is the fractional mass 
concentration of fuel in the gas mixture, m; that of the mixture leaving the reaction 
zone, and H is the calorific value of the fuel, the extinction condition corresponds 
very closely with the maximum volumetric heat release rate. Longwell calculated 
this volumetric heat release rate, obtaining B’ from experimental data on flame 
stabilisation in the wake of a bluff body. Avery and Hart‘ made alternative 
estimates of B’ by extrapolating experimental data on hydrocarbon oxidation and by 
substituting experimental values of laminar flame speed in the Zeldovich and Frank- 
Kamenetsky equation"®”. 


Now A=Uap,, and so the maximum heat release rate per unit volume of the 
primary reaction zone (q””) is equal to HU pu (Mu — Ms) 


When conditions are such that molecular transport processes are negligible 
compared with those of turbulent transport (this will be true at high values of 
Reynolds number), Ux is proportional to the velocity of the gas stream correspond- 
ing to blow-out conditions (Vo). Further v/a, which is proportional to the 
characteristic dimension of the reaction zone, is also proportional to that of 
the stabiliser (d), and hence 


. V 


The volumetric heat release rate may be obtained more directly in terms of 
laminar flame speed. For a given mixture ratio, characterised by mm, the heat 
release rate per unit area of a laminar flame is equal to Hp,S, (m:ma— my), where S, 
is the laminar adiabatic flame speed of the mixture referred to the temperature of the 
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unburned state. The thickness of the reaction zone for a laminar flame is assumed 
to be approximately equal to the flame thickness, defined“) as 


2k 
CPuSu 
where k=thermal conductivity of the gas mixture at a suitable mean 
temperature, 
and .  c=specific heat of the mixture at the same mean temperature. 


Thus, for a laminar flame, the mean volumetric heat release rate is 


H 
=K, Pu Su? (Mu — Mp) (3) 


where K, is a constant of order unity. 


Assuming that in a laminar flame the mean volumetric heat release rate is of 
the order (say 4) of the maximum given by equation (2) (for the same mixture ratio), 
q’” may be eliminated between equations (2) and (3) to give 


Vx0Pu 
K,—=* =K, S.? 


2 
or, in dimensionless form, =K, ‘ ‘ (4) 


a 


where 2=k/(cp,) and has the form of a thermal diffusivity. 


Equation (4) represents Longwell’s result rewritten in terms of Peclet numbers. 
That flame stability data might be correlated in this manner was first suggested by 
Putnam and Jensen’. 


2.2. Combustion with finite mixing rates 


Whereas Longwell assumed mixing rates within the reaction zone to be very 
fast, Spalding‘) considered the effect of finite mixing rates for various systems having 
characteristics in common with stabilised flames. For each of these hypothetical 
models, analysis showed that, when the Reynolds number is so high that molecular 
transport processes are completely outweighed by those of turbulent transport, biow- 
out phenomena should correlate in the following manner 


the symbols having the same significance as before. 


The reaction rate was, in all cases, assumed to be zero at temperatures 
corresponding to the fully-burned and unburned states, but the results were 
independent of activation energy and of the order of the reaction. 
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2.3. Dimensional analysis 


That the results obtained from the preceding analyses are independent of the 
particular model chosen is due to the limited number of dimensionless groups which 
can be formed from the relevant variables. Therefore dimensional analysis alone 
should suffice to derive the relations of Sections 2.1 and 2.2. 


The differential equation for steady-state heat transfer when heat is being 
produced by chemical reaction may be written in the form 


et? ay t” +Zp,"f(T)=0. . (6) 


where Zf(T) is a measure of the reaction rate; f(T) is a function of T alone and Z 
contains all terms relating to the initial concentrations of fuel and oxidant and their 
diffusion coefficients, which are assumed to be equal to the thermal diffusivity. This 
assumption is exact for turbulent transport alone, and nearly true for molecular 
transport. 


To find the solution for this system in terms of non-dimensional quantities, the 
method first used by Stokes and extended by Reynolds is followed (see Ref. 8). 
The solution corresponding to conditions at blow-out of the flame is being con- 
sidered here. Since there are three groups of terms in the differential equation, this 
condition can be characterised by a relation between two independent dimensionless 
quantities. The first of these is a Peclet number, derived by dividing the term 
corresponding to convective transport, in a modified form of equation (6), by 
that corresponding to conductive transport; it is Vgod/z. The second is derived by 
dividing the term accounting for chemically liberated heat by that corresponding to 
conductive transport. This is Zp,"~'d*/a, and so the equation governing extinction 
in a system of given geometrical shape is 


(Fat) (7) 


where ¢ is an unknown function. In deriving this solution, variations in the Prandtl 
and Schmidt numbers of the diffusing components have been neglected, an assump- 
tion that, in dealing with gases, is reasonably accurate. Density changes consequent 
upon reaction have also been neglected. 


The same procedure may be applied to a steady laminar flame. Since in this 
case the flame speed S, is independent of apparatus dimension, it must be possible 
to cancel d from equation (7) when the typical velocity Vgo is replaced by S,. This 
means that the function ¢ must, for this case, be a square root, and hence 


=constant=K, (say). . (8) 
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Eliminating Z between equations (7) and (8) gives, in general, 


Returning to the consideration of equation (7) and noting that at high values of 
Peclet number (this corresponds to high values of Reynolds number also) turbulent 
processes predominate, Vo must be independent of molecular transport phenomena 
(i.e. independent of k). Equation (7) can be rewritten as 
Zee) 


The only way in which the right hand side can be independent of k is if 


is linear. 


Hence, in the particular case of high Peclet numbers, equation (9) reduces to 


a 


2.4. Comments 


None of the foregoing theories provides a complete analytical solution, since 
the value of the numerical constant in each of equations (4), (5) and (10) must be 
determined from experimental data. If consistent data are used in the evaluation, 
the constants K,, K, and K,, must be equal and the equations become identical. 


Only the homogeneous reaction zone theory makes any assumption as to the 
order of the reaction and the exact form of the reaction law. Dimensional analysis, 
for example, leads to equation (10) without recourse to such assumptions. 


In none of the three approaches to the problem is any distinction made between 
what are normally designated as two- and three-dimensional stabilisers. Thus, at 
values of stream Peclet number sufficiently high for molecular transport processes 
to be negligible compared with those of turbulent transport, equation (4) will apply 
irrespective of the shape and orientation of the stabiliser (i.e. irrespective of stabiliser 
configuration). The critical value of Peclet number below which molecular trans- 
port processes are important will, however, depend upon stabiliser configuration, as 
will also the value of the numerical constant in the equation. 


3. Experimental Results 


Numerous experiments on the extinction of flames stabilised by a bluff body 
in a pre-mixed gas stream have been made and the results, usually expressed in the 
form Vo proportional to d’, have shown lack of agreement on the value of the 
exponent b. The form of the results obtained is shown by Fig. 1. Among these 
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Fig. 1. 
Stability limit curves. 


experiments are those carried out on cylindrical rods and V-gutters perpendicular 
to the gas stream (by Scurlock"’, Haddocks’, Baddour and Carr“®’, Barrere and 
Mestre“), circular cylinders with the axis parallel to the gas stream (Longwell 
et al’», Putnam and Landry“*’, discs perpendicular to the gas stream (DeZubay™’, 
Schramm'*’) and spheres (Baddour and Carr", Weir et al®”). A survey of the 
literature has been presented by Longwell"®’. 


In the light of the analyses of Section 2, it is desirable to plot these results on a 
graph of Vgod/a against S,d/z. This is done in Fig. 2. The data of as many 
investigators as possible has been included, such data being taken from the actual 
results where available, or from curves, such as Fig. 1, published in the literature. 
Since flame speed measurements for the gas mixtures used were not reported by all 
workers, suitable estimates have been made where necessary. The values of k 
and c used to form the Peclet numbers are those relevant at the mean temperature 
between the unburned and fully-burned state, while the density is that of the 
unburned gas stream. It was observed that at low Reynolds numbers the peaks of 
the stability limit curves lay off the stoichiometric line; adjustments to the air-fuel 
ratio have been made in these cases so that the peak of the stability limit curve 
corresponds with that of the flame speed curve. The shift is due to the differing rates 
at which oxygen and fuel diffuse into the reaction zone, the difference being under- 
Standably less pronounced as the stream Peclet number increases and molecular 
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Fig. 2. 


Correlation curves for flames stabilised by a bluff body in a pre-mixed gas stream. 


transport processes are gradually outweighed by those of turbulent transport. Direct 
measurements of the shifts of composition resulting from such preferential diffusion 
have been made by Williams and Shipman''”’, and the trend of the necessary correc- 
tions has followed the general experimental trend that they found. The points 
plotted include data relating to various stabiliser configurations, fuels and pressures. 
It is not to be expected that the points relating to different configurations should lie 
on a single curve, but the grouping of points and the general trend are quite 
pronounced. 


The points plotted in Fig. 2 fall close to one or other of two distinct curves. All 
points pertaining to rods and gutters perpendicular to the gas stream (i.e. the two- 
dimensional configuration) lie close to the left hand curve, while the points pertain- 
ing to discs, spheres, and rods parallel to the stream (i.e. the three-dimensional 
configuration) lie close to the right hand curve. Although data for the two- 
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dimensional configuration is lacking at high Peclet numbers, there is no reason to 
suppose that the difference in location apparent at low Peclet numbers is confined to 
this region alone. The two curves may be expected to be distinct over the whole 
range. For the three-dimensional configuration the slope of the curve, at high 
stream Peclet numbers, is 2, in agreement with the prediction of Section 2, but con- 
firmation that the left hand curve has the expected slope of 2 in fully turbulent 
flow is still lacking. At lower stream Peclet numbers the slope of the curves is of the 
order of 1-4, implying 


Vza0d alt (3) 1-4 
a a 


so that Vgo is proportional to d°*, which is the dependence found by Scurlock"? 
and others. 


A possible reason for the different slope of the curve at low stream Peclet 
numbers from that at high Peclet numbers is that, as the Peclet number (and conse- 
quently the Reynolds number) of the gas stream is decreased below a critical value 
depending upon the stabiliser configuration, the relative length of the recirculation 
zone gradually increases. This is attributable to a progressively less rapid rate of 
mixing of fresh unburned gas with the recirculation zone gases, as the relative 
importance of molecular transport processes increases. It has the effect of making 
the blow-out velocity proportional to a power of stabiliser dimension that is less 
than unity. In mathematical form, it is no longer justifiable to assume that U 
is proportional to V only. Instead, U must be replaced by a function of the form 
V (Pe,)’ (where £ is probably a positive number less than unity), which signifies that 
the velocity with which fresh gas enters the reaction zone is dependent not only upon 
the stream velocity but also upon the stream Peclet number. Equation (4) is thus 
modified and is now replaced by 


a 


and the slope of the curve is 2/(1+ 8). 


From the experimental results, the curves have a slope of approximately 1:4 
under these circumstances, and the value of f is thus 0-43. 


Qualitative experimental evidence which supports this explanation has been 
reported in the literature. Scurlock” observed a lack of similarity in the appearance 
of the recirculation zone in his experiments, the relative length of the zone changing 
with stream velocity. He also noticed that the point at which the boundary between 
the recirculation zone and the stream became fully turbulent gradually receded from 
the stabiliser as the stream Reynolds number (and hence Peclet number) decreased. 
A similar observation has been made by Williams and Shipman”. It is interesting 
to note that Scurlock’s results fall on the lower region of the curve, whereas the 
results of Longwell”), who reported a constant relative length of the recirculation 
zone, fall on the upper region of the curve, where turbulent transport processes 
predominate. 
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Fig. 3. 


Suggested form of stability curve for any given stabiliser configuration. 


At very low stream Peclet numbers a further change in the form of the curves 
might be expected, for it is probable that no flame will exist, irrespective of the rate 
at which fuel enters the reaction zone, when the stabiliser size becomes of the same 
order as the quenching distance. Since quenching distance is proportional to 2/S., 
this lower limit will occur at a constant value of S,d/a, and the curves of Fig. 2 
should be vertical at this value. A similar effect occurs in the curve of Putnam and 
Jensen”, who correlated flash-back data on a Peclet number basis. 


Figure 3 shows the general appearance of the curves correlating stability limit 
data, the dotted extension representing the suggested form at very low Peclet 
numbers when quenching distance effects predominate. 


4. The Non-Adiabatic Problem and a Method of Solution 


In the theories discussed in the previous section of this paper no account has 
been taken of heat losses from the reaction zone, except those associated with the 
mixing process (either molecular or turbulent) between burned and unburned gases. 
In these mixing processes it is a close approximation to the truth, and in addition a 
mathematical convenience, to assume that the loss of heat is so balanced by the 
gain of fresh combustible that the enthalpy remains constant. It is then possible to 
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treat the reaction rate per unit volume as a function of temperature, pressure and 
the initial state of the mixture alone. When, however, heat is lost from the zone by 
radiation, or by conduction and convection to solid boundaries, the enthalpy is no 
longer constant and the composition of the reaction zone gases is no longer uniquely 
associated with temperature. The method of calculation must therefore be modified 
so that gas composition is dissociated from temperature, and this can be done by 
following a method used for the analogous problem of solid fuel combustion"*?. 


4.1. Application of the method to the adiabatic 
reaction zone 


To illustrate the use of the method the problem of flame stabilisation in a 
homogeneous reaction zone without heat losses is first reconsidered, after which heat 
losses will be taken into account. The discussions which follow apply to any homo- 
geneous reaction zone, whether it is formed in the wake of a bluff body or in a 
combustion chamber. The reaction zone is merely a prescribed volume in which 
there is intense mixing and chemical reaction. 


Considering such a zone, there must exist a balance between the masses of fuel 
entering and leaving the zone and the mass of fuel consumed in the chemical reaction 
therein. Thus mass-balance conditions are characterised by an expression of 
the form 


Uap (my — = K,.vp" (mn, md v(T) 


where p is the static pressure and the other symbols have the same significance 
as before. 


It is necessary to satisfy simultaneously a heat balance, such that the heat 
required to raise the temperature of the unburned gas mixture to that of the reaction 
zone, together with heat losses to the surroundings, is equal to the heat liberated by 
the chemical reaction. For a lean mixture the heat-balance condition is 


. 


where H=heat of combustion of the fuel, 


Q= heat losses per unit mass of the gas mixture. 


A graphical procedure is now used to obtain the solutions of the simultaneous 
equations (12) and (13). The solutions so obtained represent the equilibrium 
conditions. 


Neglecting heat losses for the moment (i.e. @=0), equation (13), when plotted 
on rectangular co-ordinates of 7 and m;, is a straight line for any given mixture, 
characterised by m,, and initial temperature T,. This is referred to as the heat- 
balance line. Equation (12), when plotted on the same co-ordinate system, produces 
a family of curves with Ua/(vp"~') as parameter, these being referred to as mass- 
balance curves. The shape of the curves may be explained in the following manner. 
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Fig. 4. 
Heat- and mass-balance curves for the adiabatic reaction zone. 


fu 


Starting at the right hand end, if the fuel concentration is to be even slightly reduced 
by reaction, an appreciable temperature is necessary; the curve is therefore nearly 
vertical at a fuel concentration of my. However, the rapid increase of reaction rate 
with temperature ensures that thereafter quite large reductions in fuel concentration 
can be achieved with a relatively slight increase in temperature; the curve therefore 
flattens out. On the left it once again becomes steep because, as m;—> 0, there is so 
little fuel left that a very high proportion of collisions between fuel and oxygen 
molecules must be effective and hence the temperature must again increase rapidly 
if the fuel concentration is to be further reduced. 


The intersections of the mass-balance curves with the heat-balance line represent 
equilibrium conditions. Their location is shown in Fig. 4. In general the heat- 
balance line and any one mass-balance curve intersect at three points, the lowest of 
which is almost identical with the initial conditions and which, although certainly 
stable, is uninteresting since there is no flame. The highest intersection represents a 
condition of stable burning, while a little consideration will show that the inter- 
mediate intersection corresponds to a condition of unstable equilibrium. 


As the rate at which combustible gas entering unit volume of the reaction zone is 
increased, the parameter Ua/(vp"~') is similarly increased and the temperature at 
which stable burning occurs gradually falls farther and farther down the heat-balance 
line, until finally the mass-balance curve becomes tangential to this line. At this 
value of Ua/(vp"~') the stable and unstable equilibria coincide and the flame is 
on the point of blow-out. Any further increase in the rate of entry of combustible 
gas to the reaction zone lifts the mass-balance curve above the heat-balance line; 
there is no longer any intersection corresponding to stable burning and the flame is 
suddenly extinguished. The critical value of Ua/(vp"~'), above which no flame can 
exist, corresponds to the stability limit usually determined, and it is for this value of 
Ua/(vp"~*) that the volumetric heat release rate reaches a maximum. There will, 
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of course, be a separate series of curves for each set of initial conditions. The fore- 
going treatment is similar to that used by Longwell’, except that, for convenience 
in later discussion, a graphical procedure has been used to obtain the solution. 


The results indicate that there is no reason to suppose that stable burning 
cannot be maintained at indefinitely low pressures (or indefinitely high altitudes), 
provided that the rate of entry of combustible gas into unit volume of the homo- 
geneous reaction zone can be made sufficiently small, so that Ua/(vp"~') is less than 
the extinction value for the mixture in question. Experience with some combustion 
systems suggests, however, that this is not so. The following theory of the effect 
of heat losses has been developed in order to explain the discrepancy. 


4.2. Form of the solution when heat losses occur 


Firstly, heat losses by thermal radiation are considered, since these increase in 
relative importance as the pressure is reduced and are, moreover, calculable without 
exact knowledge of the shape of the reaction zone. Discussion of other losses is 
deferred until Section 6. 


The heat loss by thermal radiation from a zone of volume ~ to cold surroundings 
is proportional to pvR’E(T), where €(T) is a function increasing rapidly with 
temperature (approximately proportional to T*) and R’ is the radiation heat-loss 
factor, being dependent upon the composition and initial concentration of the 
unburned mixture and insensitive to pressure. The direct proportionality to pv is 
equivalent to assuming that all the radiation escapes from the gas without re- 
absorption, an assumption which is close to the truth for the relatively small 
dimensions and pressures under consideration. The proportionality factor R’ is so 
chosen that the heat loss per unit mass of the mixture is 


puR’€ (T) ie, UREM 
Uap ° Ua 


and equation (13) becomes 


et? ~T)+ —H (my —m)). 


Plotting equation (14), as before, gives, instead of the single heat-balance line 
shown in Fig. 4 for the adiabatic zone, a family of curves with Ua/(vR’) as 
parameter, as shown in Fig. 5. These curves droop farther and farther below the 
adiabatic line as Ua/(vR’) is decreased; there will, of course, be a separate family 
for each given set of initial conditions. The mass-balance equation (12), and con- 
sequently the mass-balance curves, remain unchanged by heat losses, since the same 
temperature is required to achieve a given amount of reaction, no matter how 
difficult it may be to attain that temperature. Stable burning conditions are now 
represented by the upper temperature intersections of the mass-balance curves with 
their corresponding partners of the heat-balance family (i.e. curves with the same 
value of Ua/v must correspond). Which pairs of curves correspond depends upon 
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Fig. 5. 
Heat- and mass-balance curves for non-adiabatic reaction zone. 


the ratio of the two parameters (i.e. p/R’) and the dashed curves in Fig. 5 
represent the loci of stable-burning intersections, each locus being valid for a single 
value of p/R’. 


It will be noticed that each locus is not unlimited in extent. It terminates at the 
upper end when the value of Ua/(vp"~') is so great that reaction can no longer 
keep pace with the rate at which fresh combustible enters the reaction zone; this is 
the type of extinction discussed in Section 4.1. It also terminates at the lower end 
when a mass-balance curve is again tangential to its heat-balance partner. Any 
further decrease in the rate at which gas enters the zone results in the zone temperature 
falling below that necessary to react the required mass of combustible; no stable- 
burning intersection exists and the flame is suddenly extinguished. This implies that 
the rate at which heat is being lost from the zone by thermal radiation is comparable 
with the rate at which heat is being liberated by chemical reaction. Further, as the 
value of p/R’ decreases, the length of the locus becomes smaller until eventually, at 
a finite value of p/R’, no stable-burning intersection exists. Thus, when thermal 
radiation losses are taken into account, there is a limiting value of p/R’ below which 
no flame can be stabilised. 


To see how the limiting values of Ua/v vary with p/R’, a cross-plot of the 
terminal values of the loci occurring in Fig. 5 gives the curve of Fig. 6, which is 
valid for the same set of initial conditions as is Fig. 5. The stable-burning region 
is that enclosed between the two arms of the curve. It is now clear how thermal 
radiation losses affect the stability limits, since for the adiabatic zone the curve of 
Ua/v against p, both plotted logarithmically, is merely a straight line with a slope of 
1/(n—1) (ie. the asymptote of the right hand arm of the curve) as shown by the 
dashed line in Fig. 6. Where the adiabatic and non-adiabatic curves are coincident it 
is justifiable to neglect radiation losses; where the two curves are different (i.e. 
at low pressures and low values of Ua/wv), this is no longer permissible. Thus, 
whereas for the adiabatic reaction zone there is no low pressure limit, for the non- 
adiabatic reaction zone there does exist a critical pressure below which no flame 
can survive. 
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Fig. 6. 


Stability limit-pressure curves for adiabatic and non-adiabatic reaction zones for given 
initial conditions. 


5. Particular Calculations 


To determine the value of this critical pressure a series of calculations were 
made for initial conditions typical of practical combustion systems; curves such as 
Fig. 7 were obtained. Fig. 8 represents the, results of cross-plotting the terminal 
values of the loci. Of course, a numerical solution can be obtained only by using 
numerical values for the terms in equations (12) and (14) and the values used are 
discussed briefly. All the calculations refer to propane gas, as a representative 
hydrocarbon, in an air mixture at an initial temperature of 289°K. 


The reaction was assumed to be the bi-molecular (i.e. n=2) collision of fuel 
and oxygen, following Zeldovich and Frank-Kamenetsky"*’, Longwell et al:?”, and 
Avery and Hart". Correspondingly the function of temperature ¥(7), in equation 
(12), was assumed to be T~*/*e-"/®, This differs from that assumed by 
Longwell) only in that he assumed the collision-factor to be independent of 
temperature, whereas the closer approximation of assuming it to be proportional 
to T! has been used here. The activation energy E was taken to be 40,000 cal./mol., 
which is the value obtained by Dugger®” from laminar flame speed data, Spalding'*” 
from the extinction of diffusion flames and Longwell”) from a combustion chamber 
in which the reaction zone was substantially homogeneous. The function of fuel 
concentrations, (mu, m;), was derived from considerations of stoichiometry. 


There is, as yet, no reliable data for the value of the collision-factor and the 
constant K,, was given a value such that, for the stoichiometric mixture at 
atmospheric pressure and the condition of blow-out, the volumetric heat release 
rate was 10° cal./(cm.* sec.), to agree with the figure calculated from laminar flame 
data for these conditions by means of equation (3). This value is about the same as 
that calculated by Avery and Hart‘ for the theoretical maximum of a homogeneous 
reaction zone operating under the same initial conditions, and it also agrees with the 
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Fig. 7. 
Heat- and mass-balance curves for a non-adiabatic reaction zone: stoichiometric mixture of 
propane-air; initial temperature 289°K. 
(The figures marked within the axes refer to pressure in atmospheres.) 


value obtained by Longwell’) in his adiabatic homogeneous reactor. The 
temperature function £(7) in equation (14) was taken as T* and a value of R’ 
appropriate to each mixture ratio considered was calculated by using the data of 
Hottel and Manglesdorf‘*’ and assuming, as mentioned earlier, that all the radiation 
escapes and that only CO, and H,O radiate. The volumetric rate of heat loss by 
thermal radiation from the products of combustion of a stoichiometric mixture at 
atmospheric pressure and flame temperature, for example, was calculated in this way 
to be approximately 0-75 cal./(cm.* sec.). 


Figure 7 shows the high temperature region of the series of heat- and mass- 
balance curves appropriate to a stoichiometric mixture for the foregoing conditions. 
while Fig. 8 is a cross-plot of the limiting conditions obtained from this and three 
similar series of curves. The temperature at which high-velocity extinction of the 
flame occurred was found to be approximately 80 per cent. of that for the fully- 
burned gas (the value varied slightly with mixture ratio and with pressure); this is in 
agreement with the value obtainable for the adiabatic zone from the equation of 
Avery and Hart’ and the experimental value obtained by Longwell’. The 
minimum stable pressure for the stoichiometric mixture was found to be of the 
order of 5 x 10-* atm., occurring at a value of Ua/v approximately equal to 2 sec.~’. 
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Fig. 8. 
Stability limit-pressure curve for non-adiabatic reaction zone : propane-air mixture; 
initial temperature 289°K. 


Further, this critical value of Ua/v turned out to be almost independent of mixture 
ratio. Fig. 8 shows curves for hydrocarbon-air mixtures of the same initial 
temperature but with mixture ratios leaner and richer than stoichiometric. It will 
be noticed that they are similar in form and horizontal position to, but are higher 
than, the stoichiometric curve. Thus mixtures with an air-fuel ratio of 28 : 1 and 
8:1 can only burn at pressures above 0-1 atm.; while very weak and very rich 
mixtures (i.e. a/f=36 : 1 approximately or a/f=5 :1 approximately) can only just 
burn even at atmospheric pressure. ’ 


6. Effect of the Assumptions on Fig. 8 


It is now necessary to examine the way in which the foregoing results depend 
upon the assumptions made in, and factors omitted from, the analysis. Owing to the 
limited amount of fundamental data at present available, only a qualitative discussion 
of the effects is possible. 


Convective heat losses from the reaction zone to the solid boundaries of the 
combustion chamber have so far been neglected. The magnitude of these losses 
depends upon the geometry of the combustion system. For the case of turbulent 
flow in the chamber, they should be expressible in the form 


K (Udp)-°? (T —Tw) 


per unit mass of entering combustible, where T, is the mean temperature of the 
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Fig. 9. 
(a) Form of curves when convective losses are also considered. 
(b) Form of curve when chemi-luminescent radiation losses are also considered. 
(c) Form of curve for reduced heat release rate or non-homogeneous reaction zone. 
(d) Form of curve for 1<n<2. 


solid boundary and d is a characteristic dimension of the combustion chamber. This 
term should be included in equation (14). The effect of such inclusion may be 
regarded as replacing the specific heat c in equation (14) by an effective specific 
heat c’, equal to c+ K (Udp)-°*. (If the flow is laminar, i.e. at very low U, c’ is 
equal to c+K(Udp)}~'.) For any given mixture ratio the curve of Fig. 8 becomes, 
for particular values of d, the dashed curves of Fig. 9(a). The curve for a particular 
d will have its minimum at a pressure corresponding to that giving a quenching 
distance* of the same order as d. Thus, the minimum pressure at which a given 
combustible mixture will burn will be increased, the increase becoming larger as the 
chamber becomes smaller, while the critical value of Ua/v at which this minimum 
occurs will be slightly decreased. 


The foregoing calculation of the radiation losses considered only thermal 
radiation of CO, and H,O, and neglected chemi-luminescent radiation, for which 
there is, as yet, insufficient data. If chemi-luminescent radiation is appreciable the 
value of R’ will be increased and the curve of stability limit against pressure will be 


*The quenching distance may be defined as the diameter of a cylindrical tube containing a 
pre-mixed combustible gas which is just too small to allow a flame to propagate along the 
tube, irrespective of the flow velocity of the gas. 
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displaced in the direction indicated in Fig. 9(b), the departure from the calculated 
curve being confined to the region where radiation losses are significant. Thus 
chemi-luminescent radiation losses (or, for that matter, an increased thermal 
radiation loss) will increase the minimum pressure at which a given mixture will 
burn and will also increase the critical value of Ua/v. 


As mentioned in Section 5, the calculations were made assuming that the 
maximum heat release rate per unit volume of the reaction zone for a stoichiometric 
mixture at atmospheric pressure is 10° cal./(cm.* sec.); it was also pointed out there 
that this is a fairly reliable estimate if the reaction zone is homogeneous. Since Ua/v 
is, for any given mixture ratio, proportional to the volumetric heat release rate, 
; the value of Ua/v at which the p=1-0 atm. line cuts the right hand branch of the 
stoichiometric curve of Fig. 8 corresponds to this heat release rate. The actual 
maximum volumetric heat release rate may, however, differ from this figure, either 
because of inaccuracies in assumptions or because of the non-homogeneity of actual 
reaction zones. In both cases the stability curve must be transposed vertically (as 
shown in Fig. 9(c)), downwards for an increase and upwards for a decrease in the 
volumetric heat release rate. This results in a change in the minimum pressure at 
which a given mixture will burn, but no change in the critical Ua/v. 


If the reaction is not of the second order, the slopes of both arms of the curves 
in Fig. 8 are changed. For an n" order reaction, the slope of the right hand arm at 
high values of Ua/v is 1/(n—1). For 1 <n< 2, which is most probable, the shape 
of the curve should be that shown in Fig. 9(d), where it is seen that the minimum 
stable pressure is decreased but the critical Ua/v remains virtually unchanged. 


7. Implications of the Results 
7.1. The altitude limit of engines 


When applied to a combustion system the parameter Ua/v has a physical 
significance, since it represents the volumetric rate of entry of combustible gas to unit 
volume of the primary reaction zone. In gas turbine engines this parameter scarcely 
varies throughout the operating range of a single engine; for the volume and entry 
area of the combustion chambers are normally fixed and, although the pressure 
varies with load, the gas velocity does not vary greatly. The operating conditions 
of an engine may therefore be represented by a vertical line on Fig. 8. 


we AV 


Combustion is stable provided that the pressure at which the chamber is 
Operating is above that corresponding to the intersection of the operating line and the 
appropriate mixture ratio curve of Fig. 8. Such intersection represents the pressure 
at which the flame is extinguished; below this pressure no flame can exist. The 
altitude corresponding to this critical pressure depends upon whether the engine is 
idling or cruising. The former is the more critical state since it gives, for a given 
altitude, the lowest combustion chamber pressure. If the altitude limit of a given 
engine is too low it may be increased by re-designing the engine, i.e. by varying the 
. value of Ua/v. Consideration of Fig. 8 shows that, if the original value of Ua/v is 

greater than 2 sec.~', it should be decreased, i.e. the chamber should be made larger; 
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TABLE I 


CALCULATED MINIMUM PRESSURES AND MAXIMUM POSSIBLE ALTITUDES FOR A 
TYPICAL ENGINE 


Limit Without Heat | WithHeat |Percentage Vari- 
Compared Losses Losses 

Ua "Pressure (atm.) 0-066 0-080 17°5 
v Altitude (ft.) 61,000 57,000 6:5 
Ua Pressure (atm.) 0-033 0:048 31:2 
v Altitude (ft.) 75,500 67,700 10-7 
Limit Pressure (atm.) zero 0-022 - 
attainable Altitude (ft.) co 85,000] 


if Ua/v is already less than 2 sec.~', however, it should be increased, i.e. the chamber 
should be made smaller. Prescription of the appropriate remedy depends, therefore, 
on correct diagnosis of the causes of extinction. 


As a practical example, consider the data for the Lucas B37 Mk. II chamber 
given by Lloyd'**). This combustion chamber operates at a value of Ua/v of approxi- 
mately 24 sec.~'. Instead of the maximum heat release rate of 10° cal. /(cm.°* sec.) used 
before for the atmospheric pressure homogeneous reaction zone, a value of 2 x 10° 
cal./(cm.* sec.) is assumed, in view of the impossibility of getting primary zone 
homogeneity without introducing a very great pressure loss; this is twice the 
maximum value reported by Mullen et al for the overall heat release of a ram-jet 
but, since the primary reaction zone only is considered here, this is probably a 
reasonable estimate. If the combustion chamber pressure at idling is 1-5 times the 
outside pressure and the minimum safe mixture ratio range is 25 per cent. above 
and below stoichiometric, the highest altitude at which the engine can idle can be 
calculated with the aid of Fig. 8. Table I shows the results of the calculations, firstly 
when heat losses are neglected and secondly when they are not. 


Suppose that, in an effort to increase the altitude limit, the size of the combustion 
chamber is doubled. The altitude limits for this larger chamber are now calculated. 
It is seen from Table I that, in the presence of heat losses, the altitude increase is not 
as great as would be expected if heat losses were neglected. The importance of heat 
losses is still more strikingly demonstrated by consideration of the maximum possible 
altitude; in the absence of heat losses there is no upper limit if the combustion 
chamber is suitably enlarged; the existence of heat losses, however, forbids operation 
of an engine of this idling pressure ratio above 85,000 ft. Published experimental 
data on altitude limits is too scanty for the predictions of Table I to be checked in 
detail. Certainly, trouble in maintaining combustion is experienced at idling above 
50,000 ft., but the lack of knowledge of how closely the actual chamber approximates 
to homogeneity makes it impossible at present to predict absolute values with an 
expectation of correctness. 
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7.2. Inflammability limits 


i A fundamental problem which awaits clarification concerns the existence of 
' limits to the range of mixture ratios which can support a propagating flame. For 
ete mixtures just outside this range it is possible to ignite the mixture by a sufficiently 
ari- }) energetic spark, but the resulting flame dies out before propagating downstream. If 

| the dimensions of the apparatus are large enough for conduction to the walls to be 
' negligible, this behaviour is inexplicable in terms of current theories of flame 
propagation, which do not recognise the existence of any “ignition temperature ” 
below which the reaction rate is zero. Accordingly a continuous range of flame 
speeds from the maximum to zero should be obtained as the mixture ratio is varied. 


The consideration of radiation losses, however, yields a possible explanation. 
- | Referring to Fig. 8, it is seen that there exists at any given pressure a mixture ratio 
for which the curve just touches the constant pressure line. If the reaction zone of a 
laminar flame can be regarded as homogeneous, this mixture is just unable to 
support a flame, because the rate at which heat is lost from the reaction zone by 
radiation is too great compared with the rate of heat liberation. Since the reaction 
zone of a flame is not homogeneous, the actual inflammability limit will lie rather 
closer to the stoichiometric mixture than is indicated by Fig. 8. For example, 
calculations indicate a lean limit for propane-air at atmospheric pressure of 1-9 per 
cent. of propane by volume; the experimental value quoted by Lewis and von Elbe®® 
is 2:12 per cent. 


8. Conclusions 


(i) Flame stabilisation phenomena will correlate satisfactorily when evaluated 
in terms of Peclet numbers. Theoretical and experimental justification for such 
correlation have been presented. 


(ii) Heat losses by thermal radiation from the primary flame region affect the 


we ' stability limits, particularly at low pressure. This leads to a minimum pressure 
be | _ below which a flame of a given combustible mixture, radiating heat to cold surround- 
tly F ings, cannot exist. This minimum occurs at a given value of the volumtetric rate of 


entry of fresh gas to unit volume of the primary reaction zone. 


on 
. (iii) The effects of other heat losses and of departures from the assumptions 
we ' made upon these critical values have been discussed qualitatively. More funda- 
ait mental data is required before this discussion can be made on a quantitative basis. 
ble Ff (iv) Possible applications of the results to engine operation and design have 
on been suggested and an explanation for the existence of inflammability limits has 
on — been put forward. 
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A Simple Approach to the Theory of 
Secondary Flows 


J. H. PRESTON, Ph.D. 
(University Engineering Laboratory, Cambridge) 


Summary: A simple method is developed for computing the trailing vorticity which 
arises when a non-uniform stream is turned. 


It is shown that, for a sudden and constant deflection of a non-uniform stream, 
no net trailing vorticity is set up in the exit flow and hence there is no secondary motion. 


In the case of an impulse cascade of finite dimensions with constant turning, it is 
found that the trailing vorticity has three distinct components—the passage vorticity 
and two components which appear as vortex sheets springing from the trailing edges 
of the aerofoils. It is shown that for small angles of deflection there is no net circulation 
associated with the trailing vorticity downstream, of the cascade, and it is inferred that 

this should still be so for large deflections. 


1. Introduction 


Secondary flows are becoming of increasing interest to the engineer. The flow 
behind wind-tunnel cascades and in the various stages of compressors and turbines 
is seriously affected by their presence. In practice these secondary flows arise when 
a boundary layer is curved in its own “ plane,” usually by an aerofoil or series of 
aerofoils abutting normal to a wall, but they also arise on swept wings, due to the 
curvature of the streamlines in the boundary layer, and on rotating aerofoils, due to 
centrifugal effects causing a bending of the streamlines in the boundary layer. Squire 
and Winter"? calculated the secondary flow behind a wind-tunnel cascade due to the 
tunnel wall boundary layer. The calculations were based on the equations of motion 
for an inviscid flow, but reasonable agreement was obtained with experiment. Haw- 
thorne’:*’ has also attacked the problem on an inviscid flow basis in a more general 
manner using vector methods. He has applied the theory to flow in pipe bends and 
to struts in boundary layers, obtaining qualitative agreement with experiment. 


The purpose of this paper is to show that the results of Squire’s and Winter’s") 
analysis and of Hawthorne’s”’ analysis can be obtained by a relatively simple 
method, making use of Kelvin’s theorem. The method is first applied to the 
simple example of a sudden deflection of a non-uniform stream and the surprising 
result is obtained that there is no net trailing vorticity downstream of the bend, and 
therefore no’ secondary flow is set up. The circular arc passage treated by Squire 
and Winter"? is next dealt with and the result derived by these authors is obtained 
much more simply. The flow through an infinite cascade of circular arc aerofoils is 
then considered. It is found for small angles of stream deflection that, at the exit 
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from the cascade, the circulation associated with the trailing vorticity contained in 
the vortex sheet springing from the aerofoil is equal and opposite to the circulation 
associated with the trailing vorticity between the blades. It is this latter which gives 
rise to the secondary flow solution found by Squire and Winter". It is suggested 
that this result should also be obtained for large stream deflections and therefore the 
form of Squire’s and Winter’s result for a circular arc passage needs modification. 
Finally the method is applied to the general case of curved non-uniform streams 
Notation 
ich 
- y measured along cascade axis perpendicular to the span 
am, Zz measured along span of cascade 
x perpendicular to yz-plane 
t is 
sty u,v,w velocity components along x, y, z, directions 
ges ¢ angle of deflection 
ion 
hat Y gap between aerofoils of cascade 
U, _ velocity at entry or exit 
q __ velocity in a passage or along a streamline 
s distance measured along a streamline 
OW 
1es n distance measured normal to the streamlines 
en 8 = constant denotes a streamline 
. a = constant denotes an equipotential line 
he 
1 1 
to h= = 
ine da/ds dB/dn 
he r radius of curvature of streamline 
on w, = dU,/dz, vorticity at entry 
W- 
al ® vorticity at exit 
id ©, component of vorticity at exit normal to streamlines 
©, component of vorticity at exit along streamlines 
aa @ angle made by vorticity vector with streamline 
le dS elementary area moving with fluid 
1e 
1g dK _ circulation around circuit enclosing dS 
id I’ circulation 
k bound vorticity 
is — 2. Assumptions 


It is assumed that the flow is steady and inviscid, although in the first instance 
viscosity will have been responsible for the production of a non-uniform stream in 
the form of a boundary layer. It must also be assumed in the case of turbulent 
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boundary layers that the Reynolds stresses can be neglected. These assumptions 
make the total pressure, 
i. 

constant along a streamline. Hence the dividing line between regions of irrotational 
and rotational flow must be a streamline. In a real flow this is obviously 
untrue, as the streamlines cross this dividing line which is the “edge” of the 
boundary layer or wake. In the calculations which follow, the “entry” flow is 
assumed to be in parallel planes and the “secondary” motion in the “exit” flow, 
which arises from vorticity components with their axes along the streamlines, is 
assumed to be negligible. This latter assumption requires the non-uniformity of the 
“entry” flow to be small and in the nature of a perturbation. In the cascade and 
curved passage problems the main flow is assumed to have fore-and-aft symmetry 
and constant turning is assumed. 


3. Theory 


The starting point of the method developed here is Kelvin’s theorem, which 
states that, for an inviscid fluid, the circulation in a circuit moving with the fluid 
remains constant provided that the density is a single-valued function of the pressure, 
that the pressure is continuous and that the extraneous forces (if any) have a 
potential. Thus the theorem will apply to a compressible fluid in the absence of 
shock waves. Use is also made of the theorem that particle paths coincide with the 
streamlines when the motion is steady, and of the Helmholtz law of vortex motion, 
which states that particles which are at any time part of a vortex filament are always 
part of the same vortex filament—a result which follows from Kelvin’s theorem. 


3.1. Sudden deflection of stream-deflector disc 


The procedure will be illustrated first for a very simple example. Imagine a 
cascade with an infinite number of closely spaced elements so designed as to deflect 
the stream of velocity U,, uniform everywhere, in planes parallel to the paper (see 
Fig. 1), through an angle «. The cascade must therefore be symmetrically placed 
relative to the incoming and outgoing streams and, in effect, it is analogous to the 
actuator disc of airscrew theory in that it produces a sudden deflection of the stream 
instead of a rise of pressure. The velocity U, is a function of the distance z, 
measured parallel to the cascade span. There is, therefore, vorticity in the oncoming 
stream of magnitude 


dU 


Wo dz 


(1) 


and the vorticity vector is at right angles to the streamlines and lies in thé xy-plane 
of the streamlines, as shown in Fig. 1. 


Consider a region of the stream of width d with a typical vortex filament at the 
position A,B, as shown in Fig. 1. Take an elementary circuit enclosing an area 
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Fig. 1. 
Sudden deflection. 


dS parallel to the streamlines and to the z-axis with its trace as shown. The 
circulation dK in this circuit is 


dK =wo,dS. : (2) 


Since the motion is steady, particle paths coincide with the streamlines and the 
circuit drifts with the fluid—the trace of dS always lying along the streamline as 
shown. The vortex filament also drifts with the fluid (since it is always composed 
of the same fluid particles) and arrives at the position A,B,. When B, moves to B, 
on the cascade, the point A, moves to A, such that 


A,A,=B,B,=d tan (¢/2). (3) 


The vortex filament is now no longer at right angles to the streamlines, but makes an 
angle # with them as shown in Fig. 1. Let B, move to B,, a distance of d tan (</2), 
then A, moves to A, through a similar distance, and 


Let the resultant vorticity be » and its components normal to and along the stream- 
lines be w, and o,. 


Now by Kelvin’s theorem w dS=w,dS. (4) 
dU 
But ®,= +, CoOt¢= —w,A,A,/(A,B,) 
= —w,2 tan (¢/2). 
Hence o,=— tan (</2). (6) 
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This then is the “secondary” trailing vorticity with its axis along the streamlines 
and to it must be added the trailing vorticity arising from the change of circulation 
about the cascade elements in the spanwise z-direction. 


The trailing vorticity is found as follows. The velocity U, may be resolved into 
x and y components 


U,=U, cos (2/2) ‘ (7) 
+V=U,sin(¢/2) fortx. . (8) 
Thus the deflector disc is a vortex sheet and about an element dy the circulation is 


The “trailing” circulation shed from an element dzdy of the vortex sheet is 


dT= z {2U, sin (2/2) dy}dz 


=2 sin (¢/2) dy dz. (10) 


Now the projection of dydz normal to the outgoing streamlines is dndz 
where dn=dy cos (¢/2). (il) 
Hence the “trailing” vorticity is 


av _ du, 
dndz dz 


Thus, adding the secondary vorticity from (6) to the “trailing” vorticity from (12), 
it is seen that the component of vorticity with its axis along the outgoing streamlines 
is zero and no secondary motion is set up downstream of the cascade. This result 
was checked by showing that it satisfied the exact equations of motion for inviscid 
flow. 


The lift dL, on an element dy dz of the deflector disc, is obtained by computing 
the change of downward momentum given to an element of mass flow pU,dn dz. 


Thus from (8), dL= pU dn dz x 2U, sin (¢/2) 
and using (11), dL= pU,? sin dy dz. (13) 
If it is assumed that aL= pU,l'dz, 


then, from (7) and (9), dL=pU,? sine dy dz, 


which agrees with (13). 
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Fig 2. 
Curved passage. 


3.2. Circular arc passage 


This problem is dealt with in the same way as in Section 3.1. 


The velocity U, is assumed to be uniform at entry and exit in planes perpendi- 
cular to the spanwise direction, and the streamlines are turned through an angle « as 
before. In the passage, following Squire and Winter“, the flow is approximated by 
a free vortex flow such that at any radius r the velocity q is given by 


when the velocity at the mean radius r, is taken to be U, (see Fig. 2). 


As in Section 3.1, a circuit is taken enclosing an area dS moving with the fluid. 
Then ©,dS is the circulation in ‘this circuit at entry, as in equation (8). The course is 
traced of a typical vortex filament A,B, through the passage (Fig. 2). Up to the 

beginning of the passage, A,B,, it remains normal to the streamlines, and on emerg- 
ing from the passage at A,B, it will be both curved and inclined to the streamlines. 
Now to pass from B, to B, a particle takes a time f,, say, given by 


on using equation (14). 


4 
| < Je 
L i= NG 
wy 
| 
ds \ 
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To pass from L to M, a time f¢ is required given by 


In time ¢, this particle has reached N, a point distant s to the right of M, given by 
ID 


In time ¢, the particle which was at P (radius r,) has reached a point R, distant s, to 
the right of Q, given by 


Hence from equations (17) and (18), 
ds’ 2er 
and =-cot¢, . . (19) 


where ¢ is the angle made by the vortex filament with the streamlines at radius r 
at the point N. At this point the velocity is now equal to U,, and so an elementary 
circuit, of area dS originally, will finally have the same area at this point. 


If, as before, © is the resultant vorticity at this point, it has components », 
normal to the streamline and —w, parallel to the streamline at this point. Hence by 
Kelvin’s theorem 


=o ,dS 


and 
dz 
Now Coto. 
Hence from equation (19), — qu, (20) 
= Ff, 
and if the passage is narrow, 
dU 


These results, giving the secondary trailing vorticity ,, are identical with those 
obtained by Squire and Winter"? using a considerably more elaborate analysis. 
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3.3. Cascade with circular arc aerofoils 


The cascade is symmetrically placed to deflect a stream that is non-uniform in 
the z-direction through an angle «< (Fig. 3). The flow between the cascade elements 
is similar to that already considered in Section 3.2 and the approximate Squire 
and Winter") solution obtained there will be assumed to hold in the passages 
between the aerofoils. As in the previous examples, the course of a typical vortex 
filament is traced through the cascade. 


VORTEX FILAMENT 
AT ENTRY 


VORTEX 
SHEETS 


Fig. 3. 


Cascade of aerofoils. 


Referring to Fig. 3, A,B, is a typical vortex filament in the entry flow, of length 
d, which lies in the xy-plane and is at right angles to the streamlines. It arrives 
at the position A,B,, where A, is the leading edge of the aerofoil and B, is d tan (</2) 
to the left of the corresponding aerofoil. The flow is assumed to meet and to leave 
the aerofoils smoothly, but nevertheless particles, just above and below the entry 
streamlines which contain the aerofoils, become separated on passing through the 
cascade, since the velocity over the upper surface of the aerofoil is higher than that 
on the lower surface. Thus, on exit from the cascade, upper and lower particles at 
B, assume positions B,’, B,, and so on, and similarly for A,, the positions A,, A,’, 
and so on; also, on exit the “A” particles are d tan (</2) in advance of the correspond- 
ing “ B” particles. Moreover, the line connecting particles, which was originally 
straight at entry, becomes curved in the region between the streamlines containing 
the aerofoils, in accordance with the calculation and assumptions made in the 
previous section. Thus, the vortex filament at exit has the appearance shown in 
Fig. 3, where between the “ aerofoil” streamlines it is curved and no longer at right 
angles to the streamlines and there is a component of vorticity along the streamlines 
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already given (approximately) by equations (20) and (20a). On the aerofoil stream- 
lines, the filament lies along the streamlines and this, by Kelvin’s theorem, must 
represent trailing vorticity infinite in magnitude—in other words trailing vorticity in 
the form of a vortex sheet springing from the aerofoil trailing edge. This vorticity is 
of opposite sign to the trailing vorticity component in the passage. In addition to this 
trailing vorticity in the form of a vortex sheet, there is trailing vorticity due to 
the spanwise variation of the “bound” vorticity at the aerofoil. This vorticity is 
also shed in the same vortex sheet springing from the aerofoil trailing edge and the 
sign is the same as the other component of this sheet. 


It is of interest to calculate the circulation contributions associated with these 
distributions of trailing vorticity for a spanwise width dz. 


The trailing vorticity in the passage between the aerofoils is given approximately 
by equation (20), i.e. 


or, averaging over the passage, by 
qu, 
dz 


The area of the passage is Y cos (</2) dz, where Y is the gap, and the circulation is 


To calculate the circulation in the trailing vortex sheet due to the vortex 
filament being bent along the streamline on passing the aerofoil, it is necessary to 
calculate the displacement of two particles—one passing over the upper surface and 
the other over the lower surface of the aerofoil. 


Let the velocities on the upper and lower surfaces of the aerofoil at any point 
distant s from the leading edge be U,+k/2 and U,—k/2, where k is the strength of 
the bound vorticity and is small compared with U,. The time taken to pass over the 
upper surface will be 


ds 
JU,+k/2 
0 
T.E 
and for the lower surface i= — (23) 
U,-k/2° 
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Therefore | ds ( ) 
1 1 


(24) 


where I" is the circulation about the aerofoil. 


Hence the distance apart of the particles is 


s=(t,-—t,)U,= 4 F (25) 
But considering a path around the aerofoil formed by the two mid-streamlines and 
two lines parallel to the cascade axis in front and behind the cascade, where the 
velocity is uniform, we obtain (see Fig. 3) 


I'=2VY =2YU, sin (</2). (26) 
Hence the distance apart of the particles is 


To calculate the strength of the trailing vortex sheet it is convenient to imagine 
the “ bound” vorticity at the aerofoils to be spread over a distance dn and that the 
paths of upper and lower particles which pass an aerofoil are separated normal to 
the streamlines by this amount. The vorticity vector now no longer lies parallel 
to the streamlines, but makes an angle ¢ with them given by 


2Ysin(</2) 
(28) 

Now by Kelvin’s theorem the component of vorticity perpendicular to the 
streamlines must be w,, and hence the component along the streamlines is 


qu, 2Y sin (</2) 


dz dn 29) 


=, 


which becomes infinite as dn—>0. The circulation I’, associated with the trailing 
vorticity is, for a spanwise width dz, 


dndz= 2Y sin(</2)dz. . (0) 
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Finally there is the circulation I’, associated with the trailing vorticity which is 
shed in the same vortex sheet due to the spanwise variation of the “ bound ” vorticity 
at the aerofoil. The circulation around the aerofoil is given by equation (26). 


Hence the trailing vortex sheet of width dz has a circulation I’, given by 


dU 
r,= 2Y sin(e/2)dz . GB 


whicii is of the same sign and equal in magnitude to I’,. Hence the net circulation 


association with the trailing vorticity in one passage and from one aerofoil in the 
form of a vortex sheet is 


+040, 


_ du, _ 

4Y sin (</2) dz dz 2<Y cos (2/2) dz 

=2 Y sin (¢/2)dz{2—ecot(e/2)} . 


It would have been expected, from Section 3.1, that this result would be true for 
all values of ¢, since any cascade of finite proportions must appear from a large 
distance to be behaving like a deflector disc, i.e. causing a sudden deflection of the 
stream. The circulation |’, associated with the trailing vorticity in the passage and 
given by equation (20) is most likely to contain the source of the discrepancy. This 
has been found by adopting the Squire-Winter"’ approximation of free vortex flow 
in the passage, and a different assumption as to the velocity distribution would have 
yielded a different form for I’,. For I’ to vanish for large deflections, the correct 
form for 1’, would appear to be 


so that the average value of the trailing vorticity at exit from the passage would be 


dU, 


w,=—4 


instead of that given by equation (20a). 


It is readily shown that the lift for unit span on each aerofoil as computed from 
the vertical rate of change of momentum of the stream is consistent with that 
given by 


L=pU,!, 
where I’ is the circulation about the aerofoil given by equation (26) and U, is the 
component of U, normal to the cascade axis. 
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3.4. General case of curved flow 
Consider a pair of adjacent streamlines (Fig. 4(a)) defined by 


= tant 
8=constan en 


dB + 8=constant 


where £ is the stream function. 


a = CONSTANT 


Fig. 4. 
General case of curved flow 


Suppose these start from a region where the velocity is uniform and equal to U,. 
Let s be the distance measured along the streamlines from a suitable datum in this 
region and let q be the velocity at the position s. Then the velocity potential is 
given by 


a= | qds. (36) 


Now a particle on the 8-streamline starting from the datum line will travel a distance 
s in a time 


On the neighbouring streamline a particle, to reach the same 2 value, will take a time 


da 


(38) 
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The difference of the times is 


@ 


(39) 


Consider now the motion of a vortex filament which at the datum position is normal 
to the streamlines. Particles which are part of the filament and which travel along 
the streamlines 8=constant and 8+d8=constant, become displaced relative to the 
normals to the streamlines 2=constant, so that the filament is no longer at right 
angles to the streamlines and the resultant vorticity has a component along the 
streamlines. 


The particle on the (8+d)-streamline is in advance of that on the £- 
streamline by 


(40) 


from equation (39), 


But 


where n is measured normal to the streamlines. Hence 


(q’) 
op” 
ain =cot¢=q 7 da. 


%9 
Hence the component of vorticity along the streamlines is given as before by 


W,=, Cot 
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Now 
where 


and applying Kelvin’s theorem, 


owing to the changing cross section of the stream tube. Hence from equation (44), 


0 


% 
dU, 
which is Squire and Winter’s”? result. 


In Fig. 4(b), where ABCD is an element dzd8 formed by a pair of streamlines 
and equipotentials and r is the radius of curvature of the f-streamline, we have 


AB=r dé=ds=hda 


—— 
r 


CD=(r+dn) do-(n+* dB) da 


(r+hdB)—— hes h ip da 


1_da_ dp 
| dn 
| 
| 
| | | | | 
_ ah 
or 
1_ 12h 
Hence, from equation (48), 
hda 1 oh 
dé ap (51) 
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Now equation (47) can be written 


and substituting from equation (51), 
Hence along a streamline hw,=-2 h? dé. (54) 
Using equation (45) qh=U, 
equation (54) becomes =-2 U, (55) 


which is Hawthorne’s”? result. 


4. Discussion 


The foregoing method of investigating the secondary flows, which may arise 
through deflecting a non-uniform stream, is seen to be very simple and involves little 
more than elementary geometry. The result that no secondary flow is set up behind 
a deflector disc, symmetrically placed, is believed to be new and the inference that 
there should be no net circulation downstream of a cascade composed of finite 
elements is of considerable importance to the wind-tunnel operator. For this to be 
true at large angles of deflection implies that the estimate by Squire and Winter"? of 
the trailing vorticity in the cascade passage at exit requires modification. It is 
thought that the discrepancy arises from the assumption by Squire and Winter of 
free vortex flow in the passage, and it would be both interesting and worth while to 
evolve a curved passage for which the velocity distribution could be expressed 
analytically so that this point could be checked. 


Although the theory is essentially a perturbation theory, there is evidence":*:°’ 
that it does yield quantitative results in fair agreement with experiment, even when 
the non-uniformity of the inflow is fairly large and secondary flows are set up in 
the outflow which are no longer small. The fact that these do not have a very 
important effect may possibly be ascribed to the theorem“? that in steady motion the 
surfaces of constant total pressure contain both the streamlines and the vortex lines. 
It is the angle between the vortex lines and the streamlines which determines the 
trailing vorticity and this angle may not be much altered by the secondary motions 
which twist the exit flow. There is, however, the danger that the secondary motions 
may transfer fluid from the boundary layer to other parts of the field, so that the 
actual distribution of vorticity downstream may be quite different from that given 
by the simple theory. 
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The neglect of viscosity and turbulent shear stresses places a limit on the useful- 
ness of the theory, particularly in its applications to problems involving pressure rise 
through the cascade. In such cases the boundary layer will thicken rapidly and may 
separate. 


Another result of the small perturbation theory and the neglect of the secondary 
motions is that in the cases dealt with it is found that the lift per unit span is 
given by 


L=pU,! 


where U, is the velocity component normal to the cascade axis and I is the circula- 
tion. This is so because the flow remains two-dimensional with these assumptions. 
In practice this relation may not be even approximately true. For instance, close to 
a wall which a cascade or aerofoil abuts, the streamlines must lie parallel to the wall. 
Hence the spanwise pressure gradients must tend to zero at the wall. Thus the lift 
will not fall to zero at the wall, but the circulation must, as the velocity is zero there. 
This probably explains the overturning produced by a cascade near a wall and the 
theory presented here cannot deal with this, as constant turning is assumed and 
indeed specified. The flow is complicated by many factors and it is important to 
specify the problem when making comparisons with experiment. For instance if 
constant turning is assumed, as in the cases dealt with, we ought to find out how to 
shape the aerofoil plan form or the twist to ensure this in the experiments, by taking 
due account of the trailing vorticity. Or we may start with a given cascade geometry 
and determine theoretically the deflection of the stream at various spanwise positions. 
Such calculations are not worth while at the present stage because the complexity of 
the problem has driven us to make some rather sweeping initial assumptions, and 
moreover the trailing vorticity in the passages and in the vortex sheets will roll up or 
mix, so that downwash calculations are of doubtful value. Nevertheless when making 
comparisons between theory and experiment the basic assumptions in the theory 
must be borne in mind. 


It is evident that further progress will be difficult and that appeal to experiment 
is necessary. Nevertheless the present theory is of undoubted value in showing 
how these secondary motions arise and in giving a qualitative picture of their 
effects*. 


5. Conclusions 


A simple method is developed for computing the secondary trailing vorticity 
which arises when a non-uniform stream is turned. The essential results derived by 
Squire and Winter" and by Hawthorne” are obtained with a minimum of analysis 
and the reason for the secondary trailing vorticity is physically evident. 


*Since the writing of this paper W. R. Hawthorne‘) has carried out a theoretical investigation 
in which the conclusions of this paper are confirmed and the analysis is extended to cover 
the asymmetrical flows. 
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It is shown that for a sudden and constant deflection of a non-uniform stream, 
no trailing vorticity is set up in the exit flow and hence there is no secondary motion. 
It is also shown that for small angles of deflection there is no net circulation down- 
stream of a cascade of finite dimensions, and it is inferred that this should be true 
also for large deflections. 
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